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We analyze the impact of a porous medium (structural disorder) on the scaling 
of the partition function of a star polymer immersed in a good solvent. We 
show that corresponding scaling exponents change if the disorder is long-range- 
correlated and calculate the exponents in the new universality class. A notable 
finding is that star and chain polymers react in gualitatively different manner 
on the presence of disorder: the corresponding scaling exponents increase for 
chains and decrease for stars. We discuss the physical consequences of this 
difference. 
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1. Introduction 

Polymer theory may serve as an archetype of an approach where the appli- 
cation of the path integral formalism leads both to a quantitative under- 
standing of a whole range of physical, chemical, and biological phenomena 
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as well as to their accurate quantitative description. 1-3 Most directly this 
is shown by the Edwards model that describes a polymer chain in terms of 
a path integral and takes into account chain connectivity and self-avoiding 
interaction. 3 A textbook derivation maps this simple two-parameter model 
to the m = de Gennes limit 4 of the 0(m) symmetric field theory. Standard 
field theoretical tools explain the origin of the scaling laws that govern poly- 
mer structural behaviour and allow to calculate the exponents that govern 
these scaling laws with high accuracy. 

One of the generalizations of the above approach extends the theory 
to describe polymers of complex structure that form networks of intercon- 
nected polymer chains. The intrinsic exponents that govern the scaling of 
a polymer network are uniquely defined by those of its constituents, star- 
like subunits known as star polymers 6 (see Fig. 1). The exponents that 



Fig. 1. Solution of star and chain polymers (blue on-line) in a good solvent immersed in 
a porous medium (black). We consider the disorder that characterizes pores (impurities) 
distribution as uncorrelated (a) or correlated (b). For the d-dimensional system, disorder 
is called long range correlated if the impurity-impurity pair correlation function decays 
for the large distances r as g(r) ~ r~ a with a < d. 

govern the scaling of star polymers are universal in that they depend on 
space dimension d and star functionality / only (the number of chains 
attached to a common center). a Currently, star polymers are synthesized 
with high functionalities and form well-defined objects with interesting in- 
dustrial, technological, and experimental applications. 7 




a. 



b. 



a Note, that for / = 1, / = 2 a chain polymer is recovered. 
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In this paper we attract attention to another recent development in the 
analysis of the scaling properties of branched polymers. 8,9 Our analysis con- 
cerns the impact of structural disorder on the scaling of polymer stars and 
chains. For lattice models, where polymers are viewed as self-avoiding walks 
(SAW) this type of disorder may be realized by forbidding the SAW to visit 
certain lattice sites, which might be interpreted as lattice dilution. For real 
polymers in solvents, structural disorder may be implemented by filling a 
porous medium by a solvent with immersed polymers as shown in Fig. 1. As 
follows from our analysis sketched below, for long-range-correlated disorder 
the polymer behaviour displays a new universality, different from that of a 
SAW on a regular (undiluted) lattice. A notable consequence is that star 
and chain polymers react in qualitatively different manner on the presence of 
disorder: the corresponding scaling exponents increase for the chains and 
decrease for the stars. We discuss to which physical consequences such a 
difference in behaviour may lead. 



2. Model and method 

The starting point of our analysis is the Edwards continuous chain model, 
generalized to describe a branched polymer structure, a star polymer. We 
describe the conformation of each arm of the star by a path r a (s), parame- 
terized by < s < S a , a = 1, 2, . . . , / (S a , the Gaussian surface of the a-th 
arm is related to the contour length of the chain) , r a (0) , corresponds to the 
central point. The partition function of the system is defined by the path 
integral: 5 

Z{S a }= j D[r 1 ,...,r / ]exp^--^J JJ 6 d (r a (0) - n(0)). (1) 

Here, the product of <5-functions ensures the star-like configuration of the 
set of / polymer chains described by the Hamiltonian: 

Q> — -L CL.O — 1 

(2) 

The first term in the r.h.s. of (2) presents chain connectivity whereas the 
second term describes an excluded volume interaction. Instead of intro- 
ducing structural disorder directly into Eq. (1), we make use of its field 
theoretical representation. The corresponding derivations are described in 
details in Refs. 5 ' 8 The relevant steps read: 

(i) we map the continuous chain model (1) onto the m = limit of 0{m) 
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symmetric field theory by a familiar Laplace transform 1 in the Gaussian 
surface variables S a to conjugated chemical potentials (mass variables) n a . 
In this procedure, the product of r5-functions in (1) is represented by a 
composite operator of a product of / m-component fields 

m 

E E ^••••M', ■■•<#• (3) 

{k} »i,...,i/=l 

Here, N ll '---> l i is a traceless tensor, (j) 1 is an i-th component of the m-vector 
field <j> and in the sum over wave vectors {k} is restricted by momentum 
conservation; 

(ii) we introduce quenched random-temperature-like disorder shifting /j, a — ► 
Ma + 5 Ha {x) by random variables S^i a (x) . These have zero mean and corre- 
lations that decay at large distances as a power law: 11 

(6n a (x)6fj, a (y)) ~ \x-y\~ a . (4) 

Here, (...) stands for the configurational average over spatially homoge- 
neous and isotropic disorder and the exponent a governs the correlation 
decay. As seen below, this leads to long range correlated disorder effects for 
a < d; 

(iii) to perform the configurational average of the free energy, we make 
use of the replica method resulting in field theoretical Lagrangean of two 
couplings uq , wq- 

C = \Y. E^ + *W + 1 E E 01 • 4) • 4) + 

a—l k a—1 {k} 

n 

^ E E i^r d (4 • c) OX • <) • (5) 

a,/3=l {fe} 

Note that the evaluation of the theory (5) involves a simultaneous polymer 
(m = 0) and replica (n — 0) limit. It is this anticipated double limit that 
allows us 8 to write the Lagrangean in terms of two coupling only. A third 
coupling appears 11 for m ^ 0. 

We apply the field theoretical renormalization group (RG) approach to 
extract the universal content of (5). In this approach, the change of the cou- 
plings u, w under renormalization defines a flow in parametric space, gov- 
erned by corresponding (3- functions f3 u (u,w), (3 w {u,w). The fixed points 
(FPs) u*,w* of this flow are the solutions to the system of equations: 
(3 u (u*,w*) = 0, f3 w (u*,w*) = 0. If a stable FP exists and is accessible, 
it determines the scaling behaviour of the polymer system. In particular, 
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for a single polymer star of / arms of equal length N in a good solvent the 
partition sum (number of possible configurations) scales as 

Z NJ ~ ffKfWf- 1 , N -> oo, (6) 

with a non-universal fugacity e M and the universal star exponent 7/. 6 The 
latter is uniquely defined by the stable FP value of the anomalous dimension 
associated with the composite operator (3). We make use of two comple- 
mentary perturbation theory expansions to calculate coordinates of the FPs 
and values of the exponents. In a first approximation we apply an expan- 
sion in 11 e = 4 — d and 8 = 4 — a which allows for a qualitative description 
of the phenomena. In a further approach we apply perturbation theory in 
the renormalized couplings u and w evaluated at fixed dimension d = 3 
for a series of fixed values of the correlation parameter a. 12 In the latter 
case we proceed within a two-loop approximation refining the analysis by 
a resummation of the divergent RG series (see Ref. 8 ~ 10 for details). The 
FP picture that arises from our calculations is shown qualitatively in Fig. 
2. Both calculation schemes display a range a\ ower (d) < a < a upper (d), of 
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Fig. 2. Fixed points G (gaussian), P (pure) and LR (long range) of the RG flow in the 
plane of the two couplings it, vi. The stable FP is shown by a square. (a) uncorrelated 
disorder a > d: P (u 7^ 0, w = 0) is stable, (b) correlated disorder a < d: LR (u ^ 0, w 7^ 
0) is stable. Crossover occurs at a = d (c.f. Fig. 1). 



values for a where the long-range-correlated FP (LR, u ^ 0, w ^ 0) is 
stable and governs polymer scaling. For a < a\ OWCI , no stable FP is found. 
This has been interpreted 8 as a collapse of a polymer coil for strongly cor- 



February 2, 2008 10:9 WSPC - Proceedings Trim Size: 9in x 6in holovatch 



6 

related disorder. For a > a uppcl the pure FP (P, u ^ 0, w = 0) is stable 
and polymer scaling is not perturbed by disorder. As far as power counting 
implies that the w-term in (5) is irrelevant in the RG sense for a > d it 
is natural to identify a up por = d. Nonetheless, in first order approximation 
of the e, ^-expansion one finds that the LR FP is stable in the unphysical 
region d < a < 2 + d/2. Our two-loop calculations at fixed d, a for d = 3 
however result in a upP er = 3 = d, ai owor — 2.2 allowing for direct physical 
interpretations . 

3. Scaling exponents 

Qualitatively, the impact of disorder can be seen already from the first order 
e, S results. Comparing the 7/ exponent (6) for the cases when the structural 
disorder is absent (or is short-range-correlated, Fig. la), 6 — 1 — e/(/ — 

3)/16, and when it is long-range-correlated (Fig. lb), 9 — 1— Sf(f— 3)/8 
we find: 

^f^J i f ) -jf = l ^(e-25), e/2<S<s. (7) 
As one can see from this estimate, the exponent difference changes sign at 

1.4 
1.2 
If 1 
0.8 
0.6 
0.4 
0.2 




12 3 4 5/ 

Fig. 3. Exponent ■yj of 3D polymer stars with / arms in a good solvent. Solid lines 
and filled squares indicate uncorrelated disorder, dashed lines S = 1.1 and open squares 
a = 2.9 correlated disorder; lines: first order e, 5-results; open, resp. filled, squares: second, 
resp. third, order resummed results for fixed d, a. 




f± = 3. The two-loop calculations slightly shift this result towards f± < 3 



February 2, 2008 10:9 



WSPC - Proceedings Trim Size: 9in x 6in holovatch 



7 

while otherwise confirming the overall picture with more accurate numeri- 
cal values for the exponents. 9 ' 10 In Fig. 3 we plot the first order e, 5-curves 
for 7/ together with the resummed two- loop estimates for d = 3. A com- 
plete account of the numerical values of the exponents is given in Ref. 9 A 
prominent feature that follows from these results is that the effect of long- 
range-correlated disorder on polymer chains and polymer stars is qualita- 
tively different. Whereas correlated disorder leads to an increase of jf for 
chain polymers (i.e. for / = 1, / = 2), the same type of disorder decreases 
the 7/ exponent for the proper star polymers (/ > 3). Below we discuss 
some possible consequences of this difference. For / > 2 we also observe 
that 7/ decreases monotonically as function of / for any valid value of a 
in the same way as it does for polymers in a good solvent: 7^ > 7^ for 
2 < /1 < f2- Furthermore, within the accuracy the data confirms a convex- 
ity from below of this function for valid values of a as well as for the e, 5 
expansion. The latter property ensures that polymer stars remain mutually 
repulsive in correlated disorder. 13 ' 14 

4. Conclusions and outlook 

It has been recognized by now that star exponents come into play for 
the qualitative description and quantitative of different phenomena, where 
statistics of branched self avoiding and random walks is involved. The exam- 
ples of such phenomena include short-range interaction of branched poly- 
mers in a good solvent), 15 diffusion-controlled reactions in the presence of 
polymers, 16 and, more generally, they are part of a multifractal description 
of diffusion limited growth in a Laplacian field. 17 Recently, star exponents 
have been used to estimate the thermal denaturation transition of DNA. 18 
Our analysis opens a way to consider the impact of long-range-correlated 
disorder on the above phenomena. 

A somewhat surprising effect that results from our calculations concerns 
the static separation in a solution of diluted chains and star polymers of 
equal molecular weight inside a porous medium. Following the estimates of 
the star exponents we predict that in a correlated medium star polymers 
will exert a higher osmotic pressure than chain polymers and in general 
higher branched star polymers will be expelled more strongly from the 
correlated porous medium. On the opposite, polymer chains will prefer a 
stronger correlated medium to a less or uncorrelated medium of the same 
density. 9 

A generalization of our approach to the case of star polymers built from 
chains of different species 14 will be presented elsewhere. 19 
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